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In this work, we investigate the characteristics of the spin-singlet state rjb of the bottomonia family 
via the radiative decays of T(nS) — > 775 + 7. The theoretical estimation of the decay widths is carried 
out in terms of the light-front quark model (LFQM). Recently CLEO and BaBar collaborations have 
measured 23(Y(3S) — > jijb) and the mass of 775. In terms of the data we fix the concerned input 
parameters in our calculations of T(nS) — > r/t + 7. A special attention is paid on the transition of 
T(5S) ->• n b + 7. The BELLE data showed that the width of T(55) -»■ T(2S, 15) + tttt is two orders 
larger than that of T(4S) — > T(2S, lS) + nn, thus some theoretical explanations have been proposed. 
Among them, it is suggested the inelastic final state interaction (IFSI) T(5S) — > BB — > Y(1S) +Tvn 
may be a natural one. If so, a similar mechanism also applies to Y(5£) — > _B(*'JjM — 5. rjb + 7, 
the precise measurement would serve as a good test whether T(5S) possess exotic components. 
Our calculation in the LFQM indicates that the rate of the direct process T(5S) — s> rjb + 7 is not 
anomalous compared to T(mS) — > r/t + 7 (m = 1,2,3,4), thus if the IFSI does apply, the rate of 
T(5S) — > Tjb + 7 should be larger than the others by orders. 
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I. INTRODUCTION 

Since the relativistic and higher-order a s corrections 
are less important for bottomonia than for any other qq 
systems, study on bottomonia may offer more direct in- 
formation about the hadron configuration and pcrturba- 
tive QCD. The spin-triplet state of bottomonia T(nS) 
and the P-states Xb(nP) were discovered decades ago, 
however its partner the singlet state rjb evaded detection 
for long time, even though much efforts were made. Many 
phenomenological researches on 775 have been done by 
some groups Different approaches result in differ- 

ent level splitting AM = T(15) - T) b (lS). In Ref. [l[ the 
authors used an improved pcrturbative QCD approach to 
get AM = 44 MeV; using the potential model suggested 
in @ Eichten and Quigg estimated AM = 87 MeV 0; in 
Ref. 3 the authors selected a non-relativistic Hamilto- 
nian with spin dependent corrections to study the spec- 
tra of heavy quarkonia and got AM=57 MeV; the lattice 
prediction is AM=51 MeV j4|, whereas the lattice result 
calculated in Ref. @ was AM = 64 ± 14MeV. Ebert 
etal. Q directly studied spectra of heavy quarkonia in 
the relativistic quark model and gave m Vb = 9.400 GeV. 

The Babar Collaboration [ic| first measured 
B(T(3S) -> jrj b ) = (4.8 ± 0.5 ± 0.6) x 10" 4 , 
M{rj b ) = 9388.91^3 + 2.7 MeV and AM = 71.4t^±2.7 
MeV in 2008 and new data were released in 2009 [Til 1 . 



More recently the CLEO Collaboration [EJ confirmed 
the observation of rjb using the database of 6 million 
T(3S') decays and assuming T{rjb) ~10 MeV, they 
obtained B(T{ZS) 7^) = (7.1 ± 1.8 ± 1.1) x 10~ 4 , 
M Vb = 9391.8 ± 6.6 ± 2.0 MeV and the hyperfine 
splitting AM = 68.5 ± 6.6 ± 2.0 MeV, whereas using 
the database with 9 million T(2S*) decays they obtained 
B(T(2S*) -> 777b) < 8.4 x 10^ 4 at 90% confidential level 
It is noted that the data of the two Collaborations 
are in accordance on M Vb , but the central values of 
i3(T(3S') — > jr/b) arc different. However, if the experi- 
mental errors are taken into account, the difference is 
still within one standard deviation. 

Some theoretical work [l3|-[l5[ is devoted to account 
the experimental results. 

In Ref. the authors studied these radiative decays 
and estimated £(T(3S) -> = 4x 10~ 4 , B(T(2S) -> 

Vb+7) = 1-5 x 10- 4 andB(T(15) 776 + 7) = 1.1 x 10~ 4 
with the mass m Vb = 9.400 GeV. Their results about m nh 
and B(T(3S) — > r/j, + 7) are close to the data. The au- 
thors of Ref. [l6[ systematically investigated the mag- 
netic dipole transition V — > P7 in the light-front quark 
model (LFQM) In the QCD-motivated effective 

Hamiltonian there are several free parameters, i.e., the 
quark mass and (3 in the wavefunction (the notation of 
/3 was given in the aforementioned literatures) which are 
fixed by the variational principle, then B(T(1S) — > 776+7) 
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1 We thank Dr. Michael Roncy for telling us the discovery history 



of rjb as the Babar collaboration published its 10-sigma discov- 
ery of the ?7i, in 2008 whereas in 2009 the CLEO collaboration 
released a re-analysis of their data and reported a 4-sigma con- 
firmation of the Babar discovery. 
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was calculated and the central value is 8.4 (or 7.7) x 10 
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FIG. 1: AM coming from different experimental measure- 
ment and theoretical work. 



It is also noted that the mass of m Vb = 9.657 (or 9.295) 
GeV presented in Ref. [l(| deviates from the data (m Vb = 
9391.8 ± 6.6 ± 2.0 MeV) and the fitted values are dif- 
ferent for singlet and triplet [l8j]. Here we are going to 
take an alternative way to fix the values of j3. 

Since experimentally, m^ b is determined by 
B(T(nS) — > rjb + 7) and a study on the radiative 
decays can offer us much information about the char- 
acteristics of rjb, one should carefully investigate the 
transition within a more reliable theoretical framework. 
That is the aim of the present work, namely we will 
use the LFQM to evaluate the hadronic matrix element 
which is governed by the non-perturbative QCD. The 
method is proven to be successful for calculating the 
transition rates of the processes where light hadrons 
exist in the final states. 

In this work, we first fix /3„'s for T(nS') in terms of their 
decay constants. Then, using the data of B(T(3S) —> 
Vb + 7) and rn Vb , we determine (3 Vb . With the parameters 
being fixed, we are able to estimate the rates B(T(1S) — > 
Vb + 7), B(T(25) -> ii b + 7), B(T(45) -> Vb + 7) and 
B(T(5S) -> Vb + l)- Since B(T(1S) -> 776 + 7) is sensitive 
to AM, the measurement of B(T(1S) r)b + 7) would 
be helpful for accurately determining the mass of 77b. 

Recently the transition rates of Y(5S) — > Y(15, 2S) + 
7T7T were measured by the BELLE Collaboration [21j and 
it was found that the widths exceed by more than two or- 
ders of magnitude the previously measured partial widths 
between lower T resonances. The authors [22| suggested 
that the re-scattering processes of T(55) -> f?(*) .£?(*) — > 



T(m5)+cr//o(980) — > T(m5)+7T7r make substantial con- 
tributions to the observable rate of the dipion transition 
of T(55) because its mass exceeds the production thresh- 
old of £?(*).£?(*), so that the intermediate bosons B^ and 
are on their mass-shell. They apply the same mech- 
anism to study the transition of T(5S) — > T(15 f ) + 77 [23[ 
and find that the re-scattering processes would enhance 
its width by almost two orders of magnitude. There 
arc different interpretations for the anomalous enhance- 
ment of the T(5iS) decays that the measured resonance 
T(10870) is a mixture of bb bound state in the 5S state 
with a hybrid bbg or a tetraquark bbqq [24| . 

Thus, we would like to further test the mechanism 
in the radiative decays of T(5S). In the radiative de- 
cay of T(5S) —> rjb + 7, the re-scattering processes 
Y(55) — > fiWfiW —5. Y(mS) + 7 also exist and one only 
needs to replace the effective vertex of B^ B^T(1S) 
and by the electromagnetic vertex B^B^j 

and B^B^rjb respectively in the diagrams given in Ref. 
[22I ]. Thus one can expect that the corresponding mech- 
anism should enhance the ratio of T(5S') — > 1% + 7- Our 
calculations show that the theoretical estimation on the 
enhancement factor strong depends on the parameter 
ffflWflWiii ( see below for more details). The future mea- 
surements on the radiative decays of T(10870) can help 
to determine if it is the T(5S') state as long as the rate of 
T(5S) — > rjb + 7 is obviously larger than that of lower res- 
onances of the family, otherwise other mechanisms may 
be more favored. Anyhow, the radiative decays would 
provide a decisive probe for the re-scattering mechanism. 

This paper is organized as follows: after the introduc- 
tion, in section II we present our calculations of the form 
factors for V — > P7 in the LFQM and the corresponding 
numerical results. In the section III we study the possi- 
ble re-scattering effects on T(55) — > rjb + 7. The section 
IV is devoted to our conclusion and discussion. 



II. T(n5) -> rj b + 7 IN THE LFQM 
A. Description of T(nS) -> rj b + 7 in the LFQM 

The Feynman diagrams describing T(nS) — > ?/b + 7 
are presented in Fig. O In this work, we calculate the 
transition rate of the radiative decays T(nS) — > rjb + 7 in 
the LFQM. 



The transition amplitude of T(nS) — > 776 + 7 can be ex- 
_ ressed in terms of tl 
is defined as 0, HH 



pressed in terms of the form factor ^ r T(nS)->i) (1 ('? 2 ) which 



(r) b (V')\J? m \r(V,h)) 



zee 



^{V ,h)q p VaJ : -T(nS)^ Vh {l )> 



(1) 



2 The different values correspond to the deferent potentials 
adopted in calculations. 



where V and V' are the four- momenta of T(nS) and rjb. 
q = V — V is the four-momentum of the emitted photon 
and e^T^/i) denotes the polarization vector of T(nS) 



3 




T(nS) 



pi 



Pi 



Vb 



-P2 
— <- 



T(nS) 




FIG. 2: Feynman diagrams depicting the radiative decay 
T(nS) -> 7 ?6 + 7. 



with helicity ft. For applying the LFQM, we first let the 
photon be virtual, i.e. leave its mass-shell q 2 = into 
the un-physical region of q 2 < 0. Then J-r( n s)^r) b (<l 2 ) 
can be obtained in the q + = frame with q 2 = q + q~ — 
q ± = — q ± < 0. Then we just analytically extrapolate 
•^ r x(nS)-i-i7b( c li) from the space-like region to the time- 
like region (q 2 > 0). By taking the limit q 2 — > 0, one 
obtains T T{nS) ^ Vb (q 2 = 0). 

By means of the light front quark model, one can ob- 
tain the expression of form factor J-YfnS)-^,,^ 2 ) [3 : 

FrinS)^^ 2 ) = e b I(mi,m 2 ,q 2 ) + e b I(m 2 ,mi,q 2 ), (2) 

where e b is the electrical charge for bottom quark , m\ = 
m 2 = m-b and 



f 1 dx f 
I(m 1 ,m 2 ,q 2 ) = J J d 2 k 



(x,k' ± )(j)(x, kj_] 
Xl M Q M' 
(kx • q±) 2 

qi 



(3) 



A = X2m\ + xim2, x = x\ and the other variables in 
Eq. ([3]) arc defined in Appendix. In the covariant light- 
front quark model the authors of Ref. [lj| obtained the 
same form factor ^ r x(?iS)-i.r; b (Q 2 )- The decay width for 
T(nS) — > rji, + 7 is easily achieved 



B. Numerical results 

In Ref. [HI the authors fixed the parameter (3 bb - = 1.145 
or 1.803 3 in the Gaussian wavefunction using the vari- 
ational method with mj = 5.2 GeV. However their pre- 
dictions on the mass m Vb = 9.657 (or 9.295) GeV and 
TO T(is) = 9.691 (or 9.558) GeV (the values in the brack- 
ets are obtained with a different potential form), which 
obviously deviate from data {m Vb = 9300 ± 20 ± 20 GeV 
and m r( i s) = 9460.30±0.26 GeV (Hi), thus we are going 
to fix /3 in an alternative way. 

First we set the b quark mass as mj = 4.64 GeV which 
was used in Ref. [l8| . Then we extract the decay constant 
h(ns) of T(nS) from the data T{T(nS) e+e") [H 
with 



T(V e+e 



47r a 2 
27 TIZ 



fv 



(5) 



In parallel, we can calculate the constant in the LFQM 
using the formula [13, EH 



fv 



4ir 3 M 



dx Id k± 



9 



^/2x(l - x)M 



-mi (mi - m 2 ) — k ± + 



mi + m 2 . 
M + mi + m 2 



(0) 



equating the two results. /3r(«s) is determined (see Table 
U). Then, we use the formula presented in Section II 
to calculate B(T(3S) — ^ ri b + 7) and compare it with 
the central value of the experimental data B(T(3S) — > 
Vb + 7) = (7.1 ± 1.8 ± 1.1) x 10-M3 or B(T(3S) -> 



% + 7 ) = (4.8 ±0.5 ±0.6) x 10" 



to fit f3 Vb and the 



corresponding values are presented in Table QT] At last 
with all the parameters we estimate B(T(1S) — > jrjb), 
B(T(2S) -> + 7), B(T(4S) -> m + 7) and B(T(55) -> 
r/b + 7) which are shown in Table [TTJ It is noted that at 
this step, we only consider the direct decay modes, but 
for T(55) — > ?/6 + 7, the re-scattering effect may play a 
dominant role as mentioned in the introduction and we 
will discuss the details in next section. 

In order to illustrate the dependence of our results on 
mfc, we re-set raj = 5.2 GeV, which was adopted by the 
authors of Ref. [l6| and fitted /3x(nS) and f3 nb again. 
Using the new parameter the B(T(nS) — > i] b + 7) is com- 
puted and the result is also listed in Table [TTJ From the 
Table HH we find 

1. The predicted branching ratios B(T(1S) —J- i] b + 7) 
and B(T(4S) — > rjb + 7) are not sensitive to m& and /3, 



r(T(nS)^ % + 7 ) 



T 2 



■>lt> 



(0),(4) 



where a is the fine-structure constant and mx( n s)i 
are the masses of T(nS) and rjb respectively. 



The authors of Ref. 1161 used the harmonic oscillator and lin- 
ear potential forms for the confinement term in their computa- 
tions, thus they obtained two different values for the parameter. 
Since later in our calculations we do not evaluate the spectra of 
the concerned hadrons, the concrete value of f) does not influence 
our numerical results. 
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TABLE I: The b quark mass and the parameters /3 in the wavefunction in the unit of GeV. 





/3r(is) 




Pr(3S) 


/?T(4S) 




4.64 


1.301 


0.924 


0.795 


0.657 


0.673 


5.20 


1.257 


0.894 


0.769 


0.635 


0.652 







TABLE II: 


The branching 


ratio of T(nS) — > 7775. 








B(T (IS) ^77^+7) B{T{2S) ->■ Vb - 


r 7) B(T(4S) ^770+7) S(T(55) ^77^+7) 


4.64 


0.915" 


2.25xl0" 4 


4.91 xl0~ 4 


4.66xl0" b 


9.27xl0" s 




0.901 6 


2.52xl0" 4 


9.21 xlO" 4 


3.86xl0 -6 


6.97xl0 -8 


5.20 


0.885" 


1.83xl0~ 4 


2.27X10" 4 


4.31x10"° 


8.57xl0" s 




0.874'' 


2.05xl0" 4 


5.31xl0 -4 


3.53xl0~ 6 


6.32xl0~ s 



"fitted from CLEO data 
'fitted from BaBar data 



but B(T(2S) -> % + 7) and B(T(5S) -> 7? 6 + 7) slightly 
change as 777,5 and /3 vary within certain ranges; 

2. Our results about S(T(15) — > i]b + 7) and 
B(Y(25) — > i]b + 7) are somehow larger than that given 
inRef. 

3. If the final state interaction is not taken into ac- 
count, the branching ratio of T(55) — > 775+7 is not 
anomalous compared to Y(mS) —> 775 + 7 (™ = 1, 2, 3, 4). 
Because the BELLE Collaboration J26| found the rate of 
T(55) — > T(1S, 25) + 7T7T is anomalously large compared 
to the similar dipion transitions between lower T reso- 



J 



nances, one has reason to doubt if such anomaly would 
appear in T(55) — > 775 + 7. 

It is noted that B(T(IS) — > 775 + 7) is sensitive to 
m rib (or AM) since the decay width is proportional to 
(AM) 3 , thus as AM is small, i.e. the masses of initial 
and daughter mesons are close to each other, any small 
changes of m Vb which has not been accurately measured 
yet [23, can lead to a remarkable difference. In Fig. Owe 
display the dependence of B(T(1S) — > 775 + 7) on AM. 
Thus the accurate measurement on B(T(IS) — > rjt, + 7) 
will be a great help to determine the mass of . 



III. POSSIBLE RE-SCATTERING EFFECTS 
INDUCING A LARGE B(T(55) ^77^+7) 



reflects the inner structures of the mesons at the effective 
vertex 



As aforementioned, the re-scattering of hadrons may 
remarkably enhance the rate of T(nS) (n > 4) — > 775 + 7. 
Similar to the re-scattering effects on in the branching 
ratios of T(55) -> T (IS) + tttt H2 and T(nS)(n > 4) -> 
T(lS)rj [23|, the transitions T(nS) — > 775+7 can occur via 
re-scattering sub-processes with the intermediate states 
being ijMijM where another B^ meson is exchanged at 
the t-channel. The corresponding diagrams are depicted 
in Fig. 21 Since the other diagrams can be obtained by 
a charge conjugation transformation of B^ the 
contribution of each diagram in Figs. 0] (a-f) should be 
multiplied by a factor 2. 

Now let us calculate the amplitudes of the re-scattering 
processes which occur at the hadron level. Following 
Refs. [22], [23| we will not account for the contribution 
from the dispersive parts of the diagrams where 
from T(ns) are off-shell, but only concern the absorptive 
parts where B^ are real particles on their mass-shells. 
The off-shell effect of the meson exchanged at t-channcl 
is compensated by a monepolc form factor which is also 



T(mi,q 



_ (A + mi) 2 - mf 
(A + rm) 2 - q 2 ' 



(7) 



where q and are the momentum and mass of the ex- 
changed meson respectively. And the cutoff is set as 
A = 600 GeV [Hlg. 

The absorptive part of the amplitude is read as 



Abs, 



Pi 



327T 2 TO T(nS) 

xd[B^B^ 



#(*) £}(*)] 



dnAi[T{nS) 
■ i] b + 7] x J 7 ^,? 2 ), 



(8) 



with i = a, b, c, d, e, /. Here, dfl and pi are the solid 
angle and linear momentum of the on-shell B*-*' in the 
rest frame of T(nS), respectively. 

The effective couplings for TBB, TB*B and TB*B* 
we adopt in this work are directly borrowed from Refs. 
|22|. l23fl , and some discussions will be made in the last 
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FIG. 3: The dependance of B(T(15) -> 7^) on AM. 
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FIG. 4: The diagrams for T(nS) -»■ S w+ S w ~ -»■ ry 6 + 
7. Other diagrams can be obtained by a charge conjugation 
transformation <=t B'*'. 



section. 



9TBB = 2.5 
9tb*b =1.4 ±0.3 
Stb-b* =2.5 ±0.4. 

Following the strategy of Ref . [23| , we obtain 

C lBB = g^BBA^BB* -BWBi), 



C 



-yB*B 



9fB'B 



m B , 

x(B*J3^-B^ B*J), 



(10) 



(11a) 



(lib) 



£ 7B , B , = g~fB*B* {—A^B*" d fxB*^ 

+A^B* u d y B*} - A^B^B*^), (11c) 
Cb-b v „ = ig B 'B nb B;d^ b B^ (lid) 

C B *B* Vb = i g -^^e^d^B* l ,B*ldpr ]b . (lie) 
tub* 

With the heavy quark spin symmetry (27j. we have 



g~fBB = gyB*B = 

g Vb B*B = g Vb B»B» 



g-tB*B*, 
— 9t(\s)bb- 



(12) 



In terms of the theoretically evaluated value of 
T(B*+ -> B+j) = 0.40 ± 0.03 kcV one can fix 
g-yB'B ~ 3.5. The coupling g Vb B* b should be at the order 
of 0(1), but may vary within a reasonable range. If we 
choose g Vb B'B = q-q b B*B* ~ 1 the contributions of the 
diagrams (a), (c), (d), (f) in Fig. |Jto B(T(5S) -> jrjb) 
are approximately a few of 10~ 7 GeV and that of the dia- 
grams (b), (e) are slightly smaller than the direct transi- 
tion. If the interference between the contributions of the 
direct decay and that through the re-scattering is con- 
structive, the total width would be about a few times 
larger, but if it is destructive, the width would be very 
suppressed. However, the coupling 3T£LS)aa ma Y be as 
large as 15 as was estimated in Refs. [22|, [23[. If we use 
this value for the coupling constants g^B* b, gr) b B'B* the 
contributions of the diagrams in Fig. [4] would enhance 
the total width by more than two orders. In principle, the 
diagram (a) in Fig. 2] can contribute to T(45) — > i] b + 7, 
but the mass of T(45) is just above the threshold of 
B + B~, so that the phase space would greatly suppress 
the contribution of the diagram. 



Ctbb = g rBB r^BB^ -Bd"Bi), (9a) 
Ctb-b = £^ e ^5„T„ 

x(B*y Bi-B*if,B*J), (9b) 
£ TB , B , = g rB , B ,(-r»B*^„B*J 

+r* i B* v d v B*J - T^B^B*^), (9c) 

^ y y ^ 

where d — d — d and the coupling constants were fixed 



IV. CONCLUSION 

By studying the radiative decay of T(nS) — > 776 + 7, 
we can learn much about the hadronic structure of 
1%. Efforts have been made to explore the spin sin- 
glet 77b, in the Data-book of 2008, i% was still omitted 
from the summary table [25| . In fact, determination 
of the mass of r\ b is made via the radiative decays of 
T(nS) —> T)b + 7 [HI, and the recent data show that 
m Vb = 9388.9±l l(stat) ± 2.7(syst) by the T(3S) data 
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and = 9394.2I^j (sta£j± 2.0(syst) by the T(25) 
data [ll| . Recently, Penin [28| reviewed the progress for 
determining the mass of rjb and indicated that the ac- 
curate theoretical prediction of m rib would be a great 
challenge. Indeed, determining the wavefunction of rjb 
would be even more challenging. In this work, we are 
not going to obtain the wavefunction or even the mass 
of rjb based on the fundamental theories, such the non- 
perturbativc QCD, but using the radiative transition to 
testify the phcnomcnologically determined wavefunction 
as long as the mass is well measured. We carefully study 
the transition rates of the radiative decays which would 
help experimentalists to extract information about m nh . 
The transition rate of T(15) — > rjb + 7 is very sensitive 
to the mass splitting AM = mx(is) — m rjb ( see the text 
above), thus an accurate measurement of the radiative 
decay may be more useful to learn the spin dependence 
of the bottominia. 

Recently our experimental colleagues have made some 
progress. CLEO and BaBar collaborator measure the 
T(3S I ) — > rjb + 7 and the mass m^ b which offer an oppor- 
tunity for us to more deeply study rjb- 

Following Ref. [l6| we systematically study T(nS) — s- 
rjb + 7 in the LFQM. We take an alternative way which 
is different from that adopted in Ref. [l(| to fix the 
parameters i.e., namely we use the data to fit /3x(„s) an d 
/3 ??b in the wavefunctions. Then using these parameters 
we estimate T(15) -> r) b + 7, T(2S*) -^rj b +j, T(4S*) -> 
rjb+j and T(55) — > rjb+j- Our result indicates T(25) — > 
rjb + 7 and Y(5<S) — > rjb + 7 are more sensitive to the 
parameters than T(15) — > rjb + 7 and T(45) — > rjb + 7. 

Since the value of T(1S) —> ?/6 + 7 is sensitive to AM, 
we show the dependance of B(T(1S) —> rjb + 7) on AM. 
We hope that our experimental colleagues will conduct 
accurate measurements in the near future to determine 
the precise value of m rib . 

The anomalous largeness of the branching ratio of 
T(55) — > T(1S, 25) + 7T7T motivates a hot surf of theoret- 
ical studies [29|. It was suggested that the re-scattering 
effects may explain the unusual large branching ratio. 
This mechanism should be tested somewhere else. In 
Ref. [23J the authors evaluated the effect induced by 
the mechanism for T(5S) —> T(1S,2S) + rj and found 
that the corresponding branching ratio is also greatly 
enhanced compared to the transition among lower res- 
onances. We suggest to further test the mechanism 
at the radiative decays where the effective electromag- 
netic vertex is relatively simple. Our result which is 
obtained in terms of the LFQM, indicates the branch- 
ing ratio of T(5S I ) — > rjb + 7 is not anomalous com- 
pared to T(mS) —> rjb + 7(771 = 1,2,3,4) as long as 
the re-scattering is not taken into account. However, 
there could be a two-order enhancement in magnitude 
for T(5S) — > %+7 which is induced by the re-scattering 
effects. Thus measurement of T(55) — > rjb + 7 would be 
an ideal probe for the re-scattering mechanism which suc- 
cessfully explains the data of T(55) — > T(1S, 2S) + nn. 
This is one of the tasks of the LHCb which will be oper- 



ating very soon. 
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Appendix 

The incoming (outgoing) meson in Fig. [5] has the mo- 
mentum P^'' = + P2 where p\^'' and P2 are the 
momenta of the off-shell quark and antiquark and 

p+ = Xl P + , p+=x 2 P+, 

Pi± = x 1 P ± + k ± , P2±=x 2 P±-k ± , 

p'+ = Xl P+, p'+ = x 2 P + , 

p'i± = XiP'± + k'±, p' 2 ± = x 2 P'± - k'± 



with x\ + x 2 = 1, where xi and k±(k' ± ) are internal vari- 
ables. M and M are defined 

xi x 2 1 

Mo = jM§-( mi -m 2 ) 2 . 



The radial wavefunctions <f> related to T(nS) are de- 
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